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/  were  conducted  for  each  of  the  range-bear ing  TMA  algorithms  and  for  the  bearings- 
only  TMA  algorithm,  which  must  be  used  whenever  range  data  is  unavailable.  The 
results  showed  that  one  particular  algorithm,  -  B'  +  (B  -  0)],  was 

signi f icantly  more  accurate  than  any  of  the  others  against  non-maneuvering 
targets  and  another  algorithm,  (2,  2  x  2;  R  -  8),  performed  best  against 
maneuvering  targets.  Both  of  these  TMA  algorithms  require  range  as  well  as 
bearing  measurement  inputs.  The  optimum  weighting  factor,  K^-p,  and  weighting 
coefficient,  a,  which  minimize  the  TMA  errors  were  also  derived. 

Least-squares  TMA  algorithms  were  derived  for  the  special  case  that  the 
target's  down-range  velocity  component  is  either  known  or  assumed  to  be  known. 

One  of  these  is  the  least-squares  linear  regression  algorithm  for  the  Ekelund 
ranging  method.  •> '  ■  ;  ‘ 
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CHAPTER  1 


INTRODUCTION  AND  OBJECTIVES 

The  development  of  passive  acoustic  systems  which  can  measure  target 
range  both  rapidly  and  continuously,  after  detection  and  tracking  have  been 
achieved,  has  produced  a  requirement  for  the  derivation  and  development  of 
Target  Motion  Analysis  (TMA)  algorithms  which  can  utilize  both  target  range 
and  bearing  data.  These  TMA  algorithms  are  used  in  Fire  Control  Systems  (FCS) 
of  submarines  and  surface  warships.  Their  primary  purpose  is  to  localize,  as 
a  function  of  time,  threat  submarine  targets  in  order  to  facilitate  weapon 
aim,  acquisition  and  target  kill,  and  to  avoid  being  counter-killed  by  target 
launched  weapons.  This  requires  TMA  algorithm  solutions  at  long  standoff 
ranges  which  are  both  timely  and  accurate  against  targets  that  may  maneuver 
frequently  in  an  unpredictable  manner.  Such  accuracy  and  timeliness  can  only 
be  achieved  with  accurate  rapid  range  and  bearing  data  as  inputs  to  the  TMA 
algorithm.  However,  since  there  are  rigorous  constraints  upon  the  range  and 
bearing  accuracies  which  can  be  achieved  with  even  the  best  platform  mounted 
sensor  systems,  it  is  important  that  only  the  most  accurate  TMA  algorthms  be 
developed  and  implemented  into  the  FCSs  of  submarines  and  surface  warships. 

An  objective  of  this  study  is  to  derive  candidate  TMA  algorithms  and 
compare  their  relative  accuracies  for  typical  encounter  geometries,  which 
include  both  maneuvering  and  non-maneuvering  targets.  A  further  objective  of 
this  study  is  to  derive,  for  each  candidate  TMA  algorithm,  the  optimum 
weighting  coefficient  for  each  range  and  bearing  measurement  versus  time  and, 
where  applicable,  the  optimum  weighting  of  the  set  of  range  measurements 
relative  to  the  set  of  bearing  measurements. 


Most  of  the  candidate  TMA  algorithms  have  previously  been  derived.  Some 
have  been  implemented  and  tested  with  field  test  data.  Although  the  bearings- 
only  TMA  algorithm  has  previously  been  error  analyzed,  none  of  those  using 
both  range  and  bearing  data  inputs  has  ever  been  adequately  error 
analyzed.^-  Therefore,  it  has  previously  been  impossible  to  quantify  and 
compare  their  accuracies. 

All  recursive  type  TMA  algorithms,  such  as  the  various  Kalman  filters, 
were  rejected  as  candidate  algorithms  since  they  are  all  inherently  less 
accurate  than  their  least-squares  linear-regression  counterparts  which  are 
non-recursive  (or  batch)  type  TMA  algorithms.  Recursive  type  TMA  algorithms 
also  tend  to  be  unstable,  whereas  the  non-recursi ve  type  TMA  algorithms  are 
inherently  stable. 

The  error  analyses  showed  that  one  of  the  TMA  algorithms  was  more 
accurate  than  any  of  the  others  for  all  of  the  non-maneuvering  target 
encounter  geometries,  and  it  was  almost  as  accurate  as  any  of  the  others  for 
those  maneuvering  target  encounter  geometries  which  were  tested.  This 
algorithm,  called  the  [KQpT  (R  -  B)  +  (B  -  0)]  TMA  algorithm,  requires  a 
reasonable  estimate  of  the  standard  deviation  range  and  bearing  measurement 
errors  in  order  to  determine  the  value  of  the  weighting  factor,  Kgpj.  The 
results,  conclusions  and  recommendations  in  this  report  should  be  of 
particular  value  to  the  developers  and  the  users  of  underwater  FCS  for 
submarines  and  surface  warships. 

^McQuitty,  J.  8.  and  Gold,  B.  A.,  Error  Analysis  of  8earings-0nly  TMA,  Naval 
Surface  Weapons  Center,  NSWC  TR  31-99,  6  May  ISB'd. 


CHAPTER  2 


DERIVATIONS  OF  TMA  ALGORITHMS 

All  TMA  algorithms  are  based  upon  the  assumption  that  the  target  moves  at 
constant  velocity  (i.e.,  constant  course  and  speed)  during  the  solution  time 
interval.  Thus,  vhe  Cartesian-  coordinate  components  of  the  target  velocity 
are  constants,  and  the  target  position  components  at  time  t  are  given  by  the 
following  two  equations  (see  Figure  2-1). 

U  +  tU  =  R  cos  B  +  X  =  Rc  +  X 

-  [R  -  E(R)][c  -  E(c)]  +  X  -  E(X)  (1) 

=  Rc  +  X  -  cE(R)  -  RE ( c )  -  E(X) 

U  +  til  -  Rc  -  X  .  -cE(R)  +  RsE(B)  -  E(X) 

W  +  till  =  R  sin  B  +  Y  =  Rs  +  Y 

-  [R  -  E(R)][s  -  E(s)]  +  Y  -  E(Y)  (2) 

=  Rs  +  Y  -  sE(R)  -  R  E(s)  -  E(Y) 

W  +  tW  -  Rs  -  Y  =  — sE ( R )  -  RcE(B)  -  E(Y) 

where 

(U,  W)  are  the  initial  target  position  components  (at  time  t  =  0) 

•  « 

(U,  W)  are  the  target  velocity  components 

R  is  the  actual  target  range  at  time  t 

R  is  the  measured  target  range  at  time  t 

E(R)  is  the  range  error  at  time  t 

B  is  the  measured  target  bearing  at  time  t 

(X,  Y)  are  the  measured  platform  position  components  at  time  t 

c  =  cos  B;  E(c)  =  — sE ( B) 

s  =  sin  B;  E(s)  =  cE(B) 


Platform  Position  at  t.  is  (X.,  Y.1 
Target  Position  at  t.  is  (U. ,  W.) 

U.  =  U  +  t.  U 
M.  =  U  +  ti  W 

Initial  Target  Position  is  (U,  W) 
Initial  Platform  Position  is  (0,  0) 


2.1  THE  (2,  2  x  2;  R  -  B)  TMA  ALGORITHM 

Two  or  more  sets  of  (R,  B,  X,  Y  and  t)  measurements  obtained  at  different 

times  can  be  used  to  generate  a  set  of  simultaneous  linear  equations  that  can 

•  • 
be  used  to  solve  for  the  four  unknowns:  U  and  U,  using  Eq.  (1),  and  W  and  W, 

using  Eq.  (2).  If  the  measurements  of  (R,  B,  X,  Y  and  t)  had  no  errors,  any 

two  sets  of  measurements,  obtained  at  different  times,  would  suffice  to  obtain 

an  error-free  solution.  However,  since  these  measurements  invariably  have 

errors,  it  is  desirable  to  try  to  minimize  their  effect  upon  the  solution 

accuracy.  A  reasonable  assumption  is  that  the  measurement  parameters  (R,  B, 

X,  Y  and  t)  are  essentially  independent  and  that  each  has  random  Gaussian- 

distributed  zero-mean  errors.  If  this  assumption  is  valid,  the  least-squares 

linear-regression  method,  invented  by  Gauss,  will  give  the  most  statistically 

•  • 

accurate  solution  of  the  four  unknowns  (U,  1),  W  and  W) .  The  first  step  in 
this  method  is  to  obtain  a  weighted  sum  of  the  squares  of  two  sets  of 
equations  like  Eqs.  (1)  and  (2),  respectively,  as  follows: 

£  ai(U  +  t  U  -  Rc  -  X)*  =  £  C-cE(R)  +  R  sE{B)  -  E(X)]J  =  E*  (3) 

1-1  1=1 

£  ojfW  +  t  W  -  Rs  -  Y)*  =  £  C-sE(R)  +  R  cE(B)  -  E(Y)]J  =  E*  (4) 

1=1  1=1 

The  next  step  is  to  obtain  the  partial  derivatives  of  Eqs.  13)  and  (4)  with 
respect  to  the  unknown  parameters  and  set  them  equal  to  zero.  This  will  give 
a  minimum  variance  equation  for  each  of  the  unknowns  as  follows: 


-  Rc  -  X)‘  = 


-  Rc  -  X).  = 


2-3 


—7  T  a.  (W+tW-Rs-Y)t=2  7  a.t.  (W  +  tW  -  Rs  -  Y).  =  0 
,u  ^  1  1  L  1  '  ' 

3  i=l  i=l 


Eqs.  (5)  and  (6)  can  be  put  into  the  following  matrix  equation: 


where 


(8) 


(9a) 


n 

B, ,  =  Z  a  =  Z  a. 

11  i=l  1 

=  I  at 
B22  =  I  at2 

Fj  =  Za  (RC  +  X) 

F2  =  Zat  (RC  +  X) 

Eqs.  (7)  and  (8)  in  matrix  form  are 


F3  =  la  (Rs  +  Y) 

84  =  Zat  (Rs  +  Y) 

Eqs.  (9a)  and  (9b)  comprise  the  (2,  2  x  2;  R  -  B)  TMA  algorithm. 

The  sum  of  the  squared  errors  in  Eqs.  (3)  and  (4),  assuming  that  E(R), 
E(B),  E ( X )  and  E ( Y )  are  independent  zero-mean  random  errors,  are  respectively 

Ex2  =  y  d.  [c2E2( R)  +  R2S2E2(B)  +  E2(X)]. 


(10) 


(11) 


Ey2  -  £  0i  [s2E2(R)  +  R2c2E2( B)  +  E2(Y)]. 
i=l 

Since  these  error  squared  terms  are  not  known,  the  best  we  can  do  is  to 
express  them  statistically  in  terms  of  their  respective  variance  as  follows: 

j2(TMA)  =  i  Z  (EX  +  Ey);  s  7x  (™A)  +  ay  (TMA) 

m-<-oo  j=l 

(12) 

=  y  a.  [a2( R)  +  R2a2(B)  +  a2(X)  +  o2(Y)]j 
1-1 

This  TMA  solution  variance,  which  is  minimized  by  the  linear  regression 

algorithm,  is  the  a  weighted  sum  of  the  down-range  variances,  Z  ao  2(R),  the 

V  2  2 

cross-range  variances,  LaR  a  (B),  and  the  platform  position  variances, 

Z  a[a2(X)  +  a2 ( Y) ] .  Normally,  the  cross-range  variance  is  much  less  than 
the  down-range  variance,  especially  at  the  longer  ranges,  since  the  cross¬ 
range  error  is  proportional  to  range  while  the  down-range  error  is 
proportional  to  range  squared.  This  particular  range-bearing  TMA  algorithm 
weights  both  of  these  errors  equally,  thus  minimizing  the  down-range  residual 
variance  at  the  expense  of  the  cross-range  variance.  The  platform  position 
variance  is,  also,  small  compared  to  the  range  variance,  which  leaves  only  the 
range  variance  as  the  dominant  term  in  Eq.  (12).  This  residual  variance  can 
be  minimized  even  further  by  weighting  each  range  measurement  variance  such 
that  each  of  them  is  equal  to  the  others.  Let  the  weighting  coefficient,  a, 
be  used  for  this  equalization.  Then, 


2-5 


a2  (TMA)  =  £  a.  [a2(R)  +  R2a2(B)  +  a2(X)  +  a2(Y)].  =  £  a.a  (R).  (13) 
1=1  1=1 

2 

Each  weighted  range  measurement  variance  term,  acr  (R),  is  equal  when  a  is 

2 

inversely  proportional  to  a  (R).  If  the  Wide  Aperture  Array  (WAA)  is  used 
to  measure  target  range,  it  can  be  shown  that  the  range  error  is 

c(R)  =  C%2°-U?  since  R  =  (14) 

l  sin^  e  A 

where 


C  *  speed  of  sound  at  the  platform 
L  =  separation  between  WAA  subarrays 
e  =  relative  target  bearing 

at  =  time  delay  difference  between  the  two  target  correlation  signals  at 
the  WAA  output 

a( at )  *  standard  deviation  error  in  the  at  measurement 
R  =  target  range 

It  can  be  shown  that  ct(t),  or  ct(at),  is  a  function  of  the  signal-to-noise  at 
the  input  to  the  correlator  processor.*  The  proposed  signal-to-noise 
weighting  of  a  due  to  a(T)  or  ct(at)  is  the  following: 


2  S/N 


The  rationale  for  this  weighting,  instead  of  (J^ )  weighting,  is  due  to  the 
presence  of  random  range  and  bearing  measurement  errors  which  are  independent 
of  signal-to-noise.  These  errors  become  the  dominant  errors  whenever  the 
signal-to-noise  dependent  errors  become  small  (i.e.,  whenever  the  signal-to- 


•Ibid. 


\  -.vV-VyWIV, 
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noise  becomes  large).  Since  the  WAA  correlator  is  a  normalized  processor,  the 
height  of  the  target  correlation  signal,  which  is  a  function  of  signal-to-noise 
at  the  input  to  the  correlator,  can  be  measured  and  used  to  calculate  the 
signal-to-noise  weighting  factor,  Eq.  (15).  Then,  combining  Eqs.  (14)  and 
(15)  gives  the  following  signal-to-noise  and  geometric  weighting  factors  which 
will  make  a  inversely  proportional  to  a  (R). 


(aS/N}i  Ug)i 


(  2  S/N  /sin^  e\2 

[rrrmji  n?"7i 


(16) 


2.2  THE  (4  x  4,  R  -  B)  TMA  ALGORITHM 


Eqs.  (1)  and  (2)  may  be  combined  Into  a  single  equation  which  can  be  used 

to  derive  the  (4  x  4,  R  -  B)  TMA  algorithm.  Multiply  Eq.  (1)  by  c  and  Eq.  (2) 

by  s  and  then  add  the  two  equations  together  to  give 

c(U  +  tU)  +  s(W  +  tW)  *  R  (cc  +  ss)  +  cX  +  sY 

=  R  +  cX  +  sY  -  E(R)  -  cE(X)  -  sE(Y) 

c(U  +  til)  +  s  (W  +  tW)  -  (R  +  cX  +  SY)  =  -E(R)  -  cE(X)  -  sE(Y)  (17) 

This  is  a  single  equation  in  the  same  four  unknowns  as  in  Eqs.  (1)  and  (2), 

0  0 

namely  (U,  U,  W  and  W).  The  least-squares  linear  regression  solution  of  a  set 
of  n  equations  of  this  type,  each  for  a  different  time,  t^ ,  can  be  derived 
in  the  same  manner  as  was  used  to  derive  Eqs.  (5)  through  (9).  The  (4x4, 
R-B)  TMA  algorithm  in  matrix  equation  form  is 


A11  a12 
A12  a22 
|  A13  A14 
i  A14  A24 


A13  A14 
A14  A24 
A33  A34 
a34  A44 


u 

°1 

• 

u 

o2 

w 

°3 

w 

04 

(18) 
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where 


A11 

•  Z.C2, 

A12  3 

Z  atC  2 

*  A13  = 

Z  a  SC  ,  A^  - 

Z  atSC 

A22 

.  Z.t2c2 

’  A24 

-Z.t 

2 sc,  A33 

-  aS  , 

ZatS2 

A44 

*  Zat2S2 

°1  * 

lac  ( R 

+  cX  + 

•  sY), 

n  -  r 

atC 

(R  +  CX  + 

sY) 

o 

II 

ZaS  (R 

+  cX  + 

'  sY), 

0,  ■  Z 

atS 

(R  +  CX  + 

SY) 

The 

residual 

variance  of 

thi s  TMA 

algorithm  as 

obtained 

mean 

errors  on  the 

righthand  side 

of 

Eq.  (17) 

is 

a2(TMA)  =  I 
m 

m 

z 

n 

I  *i. 

.  [E2(R) 

+  c 

2E2(X)  + 

S2E(Y) j 

m  -*■  °° 

j=l 

1-1 

n 

*  z 

[a 

2(R)  + 

C2a2(X) 

+  s 

2cx 2  ( Y )  ] 

i=l 

n 

*  Z 

2 

a 

(R),  , 

i  31 

since  the  variance  of  the  down-range  component  of  platform  position,  c 2<j2  (X) 
+  s2o2(Y),  is  small  compared  to  the  target  range  variance,  o(R).  Mote  the 
difference  between  this  (4  x  4,  R  -  B)  TMA  algorithm  variance  and  that  of  the 
(2,  2  x  2;  R  -  B)  TMA  algorithm,  Eq.  (13).  There  are  only  down-range  error 
terms  in  Eq.  (19),  while  Eq.  (13)  contains  equally  weighted  down-range  and 
cross-range  components  of  platform  position  and  target  position  error  terms. 
Since  the  dominant  error  term  is  the  same  for  Eq.  (19)  as  for  Eq.  (13),  the 
same  weighting  coefficient,  a,  given  by  Eq.  (16),  may  be  used  with  both  the 
(4  x  4,  R  -  3)  and  the  (2,  2  x  2;  R  -  8)  TMA  algorithms. 
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2.3  THE  8EARINGS-0NLY ,  (B  -  0),  TMA  ALGORITHM 

Eqs.  (1)  and  (2)  may  also  be  combined  into  a  single  equation  which  can  be 
used  to  derive  the  bearings-only  TMA  algorithm.  Multiply  Eq.  (1)  by  s  (*  sin  B) 
and  multiply  Eq.  (2)  by  c  (=  cos  B).  Then  subtract  Eq.  (2)  from  Eq.  (1)  to 
give  the  following  equation: 

s(U  +  tU)  -  c(W  +  tlii)  -  (sX  -  cY)  =  R  (cs  -  sc)  -  sE(X)  +  cE(Y)  (20) 

=  R[c  -  E ( c ) ]  s  -  R  [s  -  E(s)]  c  -  sE ( X )  +  xE(Y) 

=  R  E(B)  -  s£(X)  +  cE(Y) 

The  least-squares  1  inear-regression  solution  of  a  set  of  n  equations  of  this 
type,  each  for  a  different  time,  t . ,  can  be  derived  in  the  same  manner  as 
was  previously  used  to  derive  Eqs.  (9)  and  (18).  The  bearings-only  TMA 
algorithm  becomes; 


A33 

A34 

"A13 

-A14_ 

y 

V 

A34 

A44 

”A14 

-A24 

• 

0 

E2 

"A13 

-A  14 

A11 

A12 

w 

-Ai4 

"A  24 

Ai2 

*22 

w 

where  the  A's  are  the  same  as  those  used  in  Eq.  (17)  and  where 
Ei  *  Z«s(sX  -  cY),  E2  =  lots  (sX  -  cY) 

E3  =  LaC(cY  -  SX),  E4  =  late  (cY  -  SX) 

=  -L  aC  (SX  -  CY),  *  -  Z  atC  (sX  -  cY) 

The  residual  variance  of  this  bearings-only  TMA  algorithm  as  obtained 
from  the  random,  zero-mean  errors  (for  a  *  1)  on  the  right-hand  side  of 
Eq.  (20)  is 

j2(TMA)  =  £  aj  [R2a2(B)  +  S2o2(X)  +  C2a2(Y)]^ 


n  ?  ? 

=  rV(B)  , 

i  *1 
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2  2 

since  the  variance  of  the  cross-range  component  of  platform  position,  s  a  (x) 
2  2 

+  c  a  (y),  is  assumed  to  be  small  compared  to  the  target  range  variance. 

Note  that  this  (B  -  0)  TMA  algorithm  variance  has  only  cross-range  error 

terms,  as  compared  to  the  (4  x  4,  R  -  B)  TMA  algorithm  variance,  Eq.  (19), 

which  has  only  down-range  error  terms.  The  weighting  coefficient,  a,  which 

minimizes  Eq.  (22)  is  that  for  which  each  weighted  cross-range  variance  term, 

-  2  9 

a  R  7  (8)  is  equal  to  the  others.  Therefore,  a  must  be  inversely  propor- 
tional  to  R  a  (B).  The  standard  deviation  bearing  error,  cr(B),  has  the 
same  signal -to-noi se  functional  dependence  as  does  7(at),  so  that  the 
signal -to-noise  weighting  given  by  Eq.  (15)  may  be  used.  Thus,  the  weighting 
coefficient  for  the  (B  -  0)  TMA  algorithm  is 


(  2  S/N  \2 

ai  “  [j  ■+ 1  -s/ty i 


(23) 


The  (B  -  0)  TMA  algorithm  can  be  used  when  range  measurements  are 
unavailable,  since  only  target  bearing  measurements  are  required  for  a  valid 
solution.  Without  any  range  measurement  data,  the  weighting  coefficient 
reduces  to  the  signal-to-noise  weighting  factor. 


r* 


2.4  THE  [K(R  -  B)  +  (B  -  0)]  TMA  ALGORITHM 

A  weighted  range-bearing  and  bearings-only  TMA  algorithm  can  be  derived 


by  adding  a 

fraction,  K,  of 

the  (4  x  4,  R 

-  B)  TMA  algo 

(B  -  0)  TMA  algorithm,  Eq. 

(21)  as  follows 

• 

• 

r  •  • 

(KAn  +  A33) 

(KA12  +  A^) 

(K  -  1)A  13 

(K  -  DA14 

(KA12  +  A^) 

(ka22  +  a44) 

(K  -  DA14 

(k  -  da24 

m 

(K  -  1)A13 

(K  -  1)AW 

(KA33  +  Ajj ) 

(KA34  +  A12) 

L 

(K  -  1)AU 

(K  -  1)A24 

(ka34  +  a12) 

(KA44  +  A 22) 

”  -1 

U 

(KDj  +  Ej ) 

. 

u 

(kd2  +  e2) 

• 

— 

w 

(kd3  ♦  E3) 

w 

i  (K04  *  e4) 

(24) 
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The  rationale  for  this  weighted  range-bearing  and  bearings-only  TMA  algorithm 
is  to  provide  an  algorithm  which  minimizes  the  properly  weighted  sums  of  the 
down-range  variance  and  the  cross-range  variance  of  the  TMA  solution.  The 
variance  of  this  algorithm  is 


;2  2 


2v  2/ 


2%  2, 


a  (TMA)  *  V  a,  [Ka  (R)  +  RS‘(B)  +  (Kc‘  +  s‘)  a  (X)  +  (Ksc  +  c")  a  (Y)] 


This  variance  is  a  minimum  when  the  weighted  down-range  variance  sum  is  equal 
to  the  cross-range  variance  sum,  or  when 

„  _  I  a  [R2a2(B)  +  S2a2(X)  +  C2o2(Y)3 

0,,T '  f  a  c.z«)  -  eVtx)  *  TTrni 

Note  that  the  numerator  is  Eq.  (22),  the  variance  of  the  bearing-only  TMA 
solution;  and  the  denominator  is  Eq.  (19),  the  variance  of  the  range-bearing 
TMA  solution.  If  the  platform  position  variance  is  small  relative  to  the 
target  position  variance,  Eq.  (26)  reduces  to  the  following: 

„  _  I  a  R2a2(B)  HV(B) 

k°pt  ’  l  .  ,*(R>  " *  ] 

If  the  WAA  is  used  to  measure  range,  <r(R)  is  given  by  Eq.  (14),  and  K  is 
approximately  the  following: 


'OPT  ~ 


IT2 a 2 ( B )  L4  sin4 e  /l2  sin2  e\2 


a  ( B ) 
a  (At) 


If  L  =  86.6  ft,  C  ■  5,000  ft/sec,  R  ■  kyds,  <j(B)  =  degrees  and  a  (at)  *  micro¬ 
seconds,  then 

KqPJ  *  s>2?,  ?.2  (B0)  2  *  0.175  at  -li-Sly  2  (29) 

IT  a  (at,  usee)  oTa) 
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where 

R  *  kyds 

at,  ct(at)  =  u sec 
a(8)  =  degrees 

Appendix  A  gives  the  KQpj  for  0/E  or  multipath  measurements. 

2 

The  weighting  coefficient,  a,  which  minimizes  Eq.  (25),  the  a  (TMA) 

2 

for  this  algorithm,  is  obtained  by  first  minimizing  a  (TMA)  as  a  function 
of  K  by  setting  K  equal  to  KQpj  as  given  in  Eq.  (27).  This  will  make  the 
first  two  variance  terms  in  Eq.  (25)  equal  to  each  other  and  if  the  other  two 
(platform  position)  terms  are  comparatively  small,  Eq.  (25)  reduces  to 


o2  (TMA)  a  I  <*. E^opt  ct2(r)  +  r2  *  (B)^ 
=  2  I  [a  R2  a2  (B)]. 


(30) 


If  the  bearing  measurements  are  obtained  from  the  WAA,  the  basic  measurement 
is  the  time  delay,  r,  to  the  center  of  the  target  correlation  signal.  The 
relative  target  bearing  as  a  function  of  t  is 


_  _  /ct\  _  L  cos  e 

e  =  cos  (  \  or  t  = ^ 

where 

L  =  baseline  separation  of  the  WAA  subarrays 
C  =  speed  of  sound  in  the  water  at  the  array 


The  relative  target  bearing  error  is 


=  Colli 
L  sin  e 


(31) 


(32) 


If  the  platform  heading  error  is  assumed  to  be  small  compared  to  a(e),  then 
7(e)  is  approximately  equal  to  a(B),  the  true  bearing  error.  Then,  the 
weighting  coefficient  which  minimi zes  Eq.  (30)  is  inversely  proportional  to 
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the  cross-range  target  position  variance.  Thus,  the  weighting  coefficient  for 
the  [K(R  -  B)  +  (B  -  0)]  TMA  algorithm  is 

a 

Comparing  this  weighting  coefficient  with  that  for  the  (2,  2  x  2;  R  -  B)  and 
the  (4  x  A,  R  -  B)  TMA  algorithms,  Eq.  (16),  shows  that  the  geometric 
weighting  for  Eq.  (33)  is 
2 

while  for  Eq.  (16)  it  is  the  square  of  this  or 
sin  e\4 

2.5  THE  (2  x  2;  R  -  B)  TMA  ALGORITHM  FOR  AM  ASSUMED  VALUE  OF  VR 

It  is  usually  possible  to  obtain  a  more  accurate  TMA  solution  if  the 
number  of  unknowns  can  be  reduced  from  four  to  three.  For  example,  if  VR 
can  be  directly  calculated  from  the  doppler  frequency  measurement  of  a  target- 
radiated  line-component  whose  frequency  is  exactly  known,  or  if  VR  is 
assumed  to  be  known  at  any  given  time,  t^,  the  unknowns  reduce  to  U,  W  and 
VCR  where 

VCR  =  vCRft  *  X)  =  -sK  U  +  CK  W  (34) 

VR  »  VR(t  =  K)  *  Ck  U  +  Sk  W 

^  =  "Sk  VCR  +  CK  VR  (35> 

W  =  ck  VCR  +  sk  VR 
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The  TMA  algorithm  in  matrix  form  is  the  following  pair  of  equations: 


B11 

> 

B12 

B22 

• 

u 

"sk  VCR 

- 

(F1  "  ck  VR  812^ 
(F2  "  ck  VR  B22) 

"Bn 

BX2 

B12_ 

B22 

• 

M 

Ck  VCR 

(F3  '  sk  VR  B12) 
(F4  *  sk  VR  B22} 

where  the  B..  and  the  are  the  same  as  those  in  Eq.  (9).  This  TMA 
algorithm  solves  for  four  unknowns,  U,  W,  and  two  values  of  V^.  It  can  be 
reduced  to  only  three  unknowns  by  solving  for  the  weighted  mean  of  the  two 
V^R  solutions. 

VCR  =  sk  (sk  VCR^  +  ck  (ck  VCR*  =  VCR  (t  =  V  {37) 


The  value  of  VR  at  time  t^ 
range,  R,  at  time  tk  is 

Rk  =  (ck  U  +  sk  W)  +  t 

*  Rj  +  t 


is  assumed  to  be  known,  and  the  solution  for 

k  VR  "  ^ck  XK  +  sk  V 
k  VR  "  \  RPLATF0RM 


(38) 


As  previously  mentioned,  the  range  rate  can  be  obtained  by  measuring  the 
difference  between  the  doppler-modified  received  frequency  and  the  known 
radiated  frequency  as  follows: 

\  *  (VR  '  RPLATFORM^K  *  (VR  '  Vp  C0S  e)K  =  "3  T“ 'KHz’  knotS  (39) 


where 

Vp  3  platform  speed 
e  =  relative  target  bearing 
Solving  this  equation  for  VR  g.,-ves 

VR  *  (R  ^  Vp  cos  e)K  (40) 
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This  algorithm  can  also  be  used  when  the  target  is  in  the  convergence 
zone  and  is  kept  there  all  during  the  solution  time  by  appropriate  platform 
maneuvers.  For  this  case,  VR  at  time  t^  is  assumed  to  be 


VR  * 


ck  XK  +  sk  YK 


(41) 


The  solution  for  range,  R,  at  time,  t^,  from  Eq.  (38),  reduces  to 

«K  =  Rj  =  ck  U  +  S]<  w  (42) 

Since  the  range  is  held  almost  constant  by  appropriate  platform  maneuvers,  the 
accuracy  of  this  range  solution  may  be  improved  by  solving  for  the  median 
range.  Whenever  the  median  range  measurement  is  used,  the  cross-range 
velocity  solution,  Eq.  (37),  should  be  modified  as  follows: 


,VCR* 


MED 


MEDIAN  (R) 


'CR 


(43) 


The  rationale  for  using  the  median  range  instead  of  the  average  range  is  due 
to  the  skewness  of  the  range  error  distribution  which  can  result  in  a  large 
average  range  bias  error,  especially  at  the  longer  ranges,  whereas  the  median 
range  is  statistically  unbiased. 


2.6  THE  (B-0)  TMA  ALGORITHM  FOR  AN  ASSUMED  VALUE  OF  VR 

The  bearings-only  TMA  set  of  equations  in  Cartesian  coordinates  is  the 
following: 
n 

£  Csi  (U  +  ti  U  -  Xj)  -  cf  (H  +  tj  W  -  =  0]  (44) 

1-1 
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Changing  the  velocity  components  from  Cartesian  to  polar  geometry,  using  Eq. 


(35),  and  assuming  that  VR  is  known,  gives 

Y.  'Si  (U  "  ti  Sk  VCR5  ‘  Ci  fW  +  ti  Ck  VCR5  *  fsi  Xi  '  Ci  Yi) 

i  =  l 

“VR  ti  (si  ck  '  ci  sk^ 

The  TMA  algorithm  in  matrix  form  is 


(45) 


A33  A34  "A13  ~A14 

U 

^E1  "  VR  fcK  A34  ‘  SK  A14^ 

A34  A44  "A14  "A24 

“sk  VCR 

^E2  -  VR  '  CK  A44  ‘  SK  A24^ 

-A13  -A14  An  A12 

• 

W 

tE3  +  VR  (cK  A14  ‘  SK  A12}- 

"A14  ~A24  A12  A22 

CK  VCR 

[E4  +  Vr  (ck  A24  -  sK  A22)] 

where  the  A„  and  the  E.  are  the  same  as  those  in  Eq.  (21).  This 
algorithm  also  solves  for  four  unknowns,  which  can  be  reduced  to  only  three 
unknowns  by  using  Eq.  (37)  to  solve  for  the  weighted  mean  value  of  the  two 
VCR  solutions.  The  solution  (U,  W  and  VCR}  plus  the  assumed  known  value 
of  VR  can  be  converted  back  to  Cartesian  velocity  coordinates  by  use  of 
Eq.  (35). 


This  bearings-only  TMA  solution,  Eq.  (46),  may  be  particularly  useful  at 
long  ranges  since  the  "range  collapse"  phenomenon  can  be  prevented  by  a 
measured  or  judiciously  assumed  value  for  VR.  The  bearings-only  TMA 
algorithm  minimizes  only  the  cross-range  errors.  This  can  cause  large  down- 
range  errors  since  they  are  not  minimized.  By  allowing  VR  to  be  set  at  an 
independently  measured  value,  or  at  some  value  that  is  apriori  known  to  be 
reasonable,  a  more  accurate  solution  will  usually  result  at  the  longer  ranges. 


TT 


■U  ■!  T 
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When  it  is  assumed  that  VR  =  0,  this  TMA  algorithm  reduces  to  a  ''east- 

squares  linear  regression  formulation  of  the  Ekeiund  bearings-only  range  and 

cross-range  velocity  algorithm.  Note  that  when  VR  =  0,  Eq.  (46)  reduces  to 

•  • 

Eq.  (21),  except  that  the  velocity  components  U  and  W  are  changed  to  -s^  V^R 

•  • 

and  c^  V^R  respectively.  Then  using  Eq.  (37),  Vj.R  =  -s^  l)  +  ck  W,  since 
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SIMULATION  ANALYSIS  OF  THE  TMA  ALGORITHMS 

A  simulation  analysis  was  conducted  using  the  TMA  algorithms  described  by 
Eq.  (9),  the  (2,  2  x  2;  R  -  B)  algorithm;  Eq.  (18),  the  (4  x  4;  R  -  B) 
algorithm;  Eq.  (21),  the  (B  -  0)  algorithm;  and  Eq.  (24),  the  [K(R  -  B)  + 

(B  -  0)]  algorithm.  The  objectives  of  this  simulation  analysis  were  to: 

•  compare  the  TMA  solution  accuracies  of  these  four  algorithms, 

•  confirm  the  validity  of  Eq.  (27),  which  derives  the  weighting 
coefficient,  K,  that  minimizes  the  variances  of  the  [K(R  -  B)  + 

(B  -  0)]  TMA  solution  errors, 

•  confirm  the  validity  of  Eqs.  (16)  and  (33)  which  give  the  weighting 
coefficients  a.,  that  minimize  the  TMA  solution  variances  of  the 
(2,  2  x  2,  R  -  B)  and  the  [K(R  -  B)  +  (B  -  0)1  algorithms, 
respectively,  and 

t  test  the  effect  of  a  target  maneuver  upon  the  TMA  solution  accuracies 
of  these  four  algorithms. 

A  number  of  statistical  methods  could  have  been  used  to  perform  the 
statistical  analysis  of  the  TMA  algorithms.  The  method  that  was  used  was 
selected  because  it  provides  a  zero-mean,  symmetrical  error  distribution  of 
the  target  data  (i.e.,  the  target  bearing  and  the  delta-tau,  at,  used  to 
calculate  target  range)  as  a  function  of  time.  In  this  method,  the  error  for 
each  bearing  is  either  plus  or  minus  a  fixed  value,  and  the  error  for  each 
delta-tau  is  also  either  plus  or  minus  a  fixed  value,  independent  of  the 
bearing  error.  A  total  of  either  four  or  six  sets  of  data  as  a  function  of 
time  was  used.  A  TMA  solution  was  obtained  for  each  and  every  possible 
permutation  of  bearing  and  delta-tau  errors  which  gave  a  total  of  24*24  =  256 


solutions  for  the  four  sets  of  data  and  2^*2^  =  4096  solutions  for  the  six 
sets  of  data.  Each  TMA  solution,  consisting  of  the  X  and  Y  components  of  both 
initial  target  position  and  velocity,  was  used  to  solve  for  target  range, 
down-range  velocity  and  cross-range  velocity  for  the  end  of  the  solution 
time.  The  mean  and  standard  deviation  of  these  256  or  4096  solutions  were 
then  calculated  and  used  to  specify  the  standard  deviation  and  bias  of  the 
target  range,  and  of  the  down-range  and  cross-range  components  of  target 
velocity  at  the  end  of  the  solution  time. 

3.1  TMA  ALGORITHMS  ACCURACY  COMPARISONS 

The  first  simulation  analysis  was  used  to  compare  the  accuracies  of  the 
(2,  2  x  2;  R  -  B)  and  the  (4x4,  R  -  B)  algorithms  for  a  selected  group  of 
typical  TMA  run  geometries  as  shown  in  Figure  3-1.  As  expected,  the  (2,  2  x  2 
R  -  B)  algorithm  invariably,  gave  more  accurate  TMA  solutions  than  the  (4  x  4, 

R  -  B)  algorithm,  which  had  horrendous  cross-range  velocity  and  TMA  bearing 
errors,  especially  at  long  ranges.  The  cause  of  these  large  cross-range 
target  position  and  velocity  errors  can  be  inferred  from  an  inspection  of  Eq. 
(18),  the  residual  variance  of  the  (4  x  4,  R  -  B)  TMA  algorithm.  Note  that 
the  only  errors  in  this  equation  are  the  down-range  errors.  The  linear 
regression  algorithm  minimizes  these  errors,  but  it  does  so  at  the  expense  of 
the  cross-range  errors,  which  can  increase  virtually  without  bounds.  For 
example,  the  cross-range  target  velocity  error  can  be  as  much  as  hundreds  of 
knots  for  reasonable  values  of  range,  range  measurement  error  and  bearing 
measurement  error.  In  addition  to  having  cross-range  errors  that  were 
hundreds  of  times  larger  those  of  the  (2,  2  x  2;  R  -  B)  algorithm,  the  (4  x  4, 
R  -  B)  algorithm  also  had  consistently  larger  range  and  down-range  velocity 
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errors.  This  results  from  the  fact  that  the  14  x  4,  R  •  B'  algorithm  must 
simultaneously  solve  for  four  unknowns,  whereas  the  (2,  2x2;  R  -  B) 
algorithm,  which  consists  of  two  sets  of  equations,  must  simultaneously  solve 
for  only  two  unknowns.  This  analysis  eliminated  the  (4  x  4,  R  -  B)  algorithm 
as  a  candidate  for  TMA. 

3.2  WEIGHTING  COEFFICIENT,  KQpT 

Next,  the  TMA  solution  accuracy,  as  a  function  of  the  weighting  factor 
constant,  K,  was  obtained  for  the  [K(R  -  B)  +  (B  -  0)]  algorithm,  Eq.  (24), 
for  values  of  K  between  0  and  1  inclusive.  Note  that  when  K  =0,  this 
algorithm  becomes  the  (B  -  0)  algorithm,  Eq.  (21),  and  when  K  =  1,  it  becomes 
the  (2,  2  x  2,  R  -  B)  algorithm,  Eq.  (9).  The  value  of  K  which  gives  the 
minimum  TMA  solution  error  is  K0pT  as  given  by  Eq.  (27).  The  typical  TMA 
run  geometries  shown  in  Figure  3-1  were  used  to  determine  the  TMA  solution 
errors  as  a  function  of  K.  The  results  showed  that  ail  of  the  TMA  solution 
errors  were  at,  or  very  near,  their  minimum  values  for  K  *  KQpj.  They  a1 so 
showed  that  for  values  of  K  between  about  0,4  KqPT  and  2.5  Kgpj,  the  TMA 
solution  errors  were  very  near  their  minimum  values.  Figure  3-2  shows  o(R)  as 
a  function  of  K.  Figure  3-3  shows  o(R),  ?(VR)  and  o(V(;R)  as  a  function  of 
K  for  the  three  most  significant  values  of  K:  K  =  0,  the  (B  -  0)  TMA 
algorithm;  K  =  KqPT>  the  TKqp-j-  (R  -  B)  +  fB  -  0)]  TMA  algorithm;  and  K  = 

1,  the  (2,  2  x  2;  R  -  B)  TMA  algorithm.  The  TMA  solution  errors  o(R),  a(VR) 

and  ° ( VqR)  are  the  standard  deviation  errors  of  the  TMA  solutions  for  R, 

VR  and  V^R ,  respectively,  at  the  end  of  the  solution  time  interval  ( i . e. , 
at  the  end  of  the  second  leg  of  platform  motion).  Note  that  the  (8  -  0)  TMA 

errors,  (K  =0),  are  less  than  the  (2,  2  x  2;  R  -  8)  TMA  errors,  rK  =  1),  for 


some  of  the  runs  and  greater  for  the  other  runs,  but  the  [Kgpy  (R  -  B)  + 

(B  -  0)]  TMA  errors  (K  =  Kgpy)  are  less  than  either  for  all  of  the  runs. 

It  is  interesting  to  compare  in  Figure  3-2  the  u(R)  versus  K  curves  of 
runs  Nos.  1  and  2.  In  run  No.  1,  the  platform  conducts  lag-lead  legs  of 
motion  and  in  run  No.  2  it  conducts  lead-lag  legs  of  motion,  against  a  target 
moving  with  the  same  cross-range  and  down-range  velocity  components.  The  (B  - 
0)  TMA  algorithm,  (K  =  0),  gives  a  ar(R)  of  11  kyds  for  run  No.  1  and  a  cr(R)  of 
only  2.35  kyds  for  run  No.  2.  This  large  percentage  difference,  (468 
percent),  in  the  a(R)  of  the  lag-lead  versus  the  lead-lag  platform  runs  for 
bearings-only  TMA  was  initially  described  in  a  previous  technical  report.* 

The  (2,  2  x  2,  R  -  B)  TMA  algorithm  gives  a  c(R)  of  3.8  kyds  for  both  runs  No. 
1  and  No.  2,  thus  indicating  that  this  algorithm  gives  the  same  TMA  accuracy, 
independent  of  whether  the  platform  motion  is  lag-lead  or  lead-lag.  However, 
the  most  accurate  TMA  solution  is  obtained  by  using  the  [(K0pT(R  -  B)  +  (B  - 
0)]  algorithm  and  conducting  lead-lag  legs  of  platform  motion,  as  in  run  No. 

2.  This  gives  a  o(R)  of  2.0  kyds.  If  the  lag-lead  legs  of  platform  motion, 
as  in  run  No.  1,  were  conducted  using  this  algorithm,  a(R)  would  have  been  3.5 
kyds,  which  is  about  75  percent  greater  than  for  lead-lag  legs.  Therefore, 
even  though  this  percentage  difference  is  considerably  smaller  than  that  for 
the  (8-0)  algorithm,  it  is  still  large  enough  to  justify  using  lead-lag  legs 
rather  than  lag-lead  legs  of  two-leg  platform  motion.  Appendix  B  shows  why 
the  lead-lag  platform  motion  tactic  gives  a  more  accurate  TMA  range  solution 
than  the  lag-lead  platform  motion  tactic. 


1 1  b  i  d . 


3.3  GEOMETRIC  WEIGHTING  COEFFICIENT,  ag 

The  weighting  coefficient,  <*,  for  weighting  each  linear  equation  relative 
to  the  others,  in  the  set  of  equations  which  comprises  the  TMA  algorithm,  is 
the  product  of  two  different  weighting  factors,  the  signal-to-noise  weighting 
factor,  Eq.  (15),  and  the  geometric  weighting  factor,  ag.  Eq.  (33), 

which  gives  the  optimum  weighting  coefficient,  a,  for  the  [KqPj  (R  -  B)  + 

(B  -  0)]  TMA  algorithm,  shows  a  geometric  weighting  factor  of 

/  .  \  7 

sin  e\  . 


The  optimum  geometric  weighting  factor  for  the  (2,  2  x  2;  R  -  B)  TMA  algorithm 
is  shown  in  Eq.  (16)  to  be 


sin  e 


or  the  square  of  that  for  the  [K(R  -  8)  +  (B  -  0)]  algorithm. 

The  five  run  geometries  shown  in  Figure  3-4  were  simulation  analyzed  in 
order  to  determine  the  validity  of  these  geometric  weighting  factors.  Six 
sets  of  data  were  used  to  obtain  each  solution,  which  gave  a  total  of  4,096 
solutions  per  run  that  were  used  to  obtain  the  standard  deviation  and  bias  of 
the  TMA  range  solution. 

The  range  bias  error  is  due  to  the  non-linearity  of  the  range  error  as  a 
function  of  bearing  error.  If  the  bearing  error  is  assumed  to  be  a  zero-mean 
random  error  with  a  symmetrical  distribution,  a  reasonable  assumption,  the 
range  error  distribution  is  skewed  (i.e.,  it  is  both  biased  and  non- 
symmetrical ) .  A  simple  example  will  show  this  non-linearity.  Let  the  range 
to  a  stationary  target  be  determined  by  measuring  the  target  bearing  from  two 
different  locations,  A  and  B,  that  are  separated  by  a  distance  L  which  is 


•Ibid. 
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Figure  3-4.  Run  Geometries  Used  to  Analyze  the  Geometric  Weighting  Factor 


assumed  to  be  small  compared  to  the  range.  For  this  case  the  actual  range  is 


R~ 


sin  ?  (eA  +  «B) 
9B  '  eA 


L  sin  e 
ae 


(46) 


If  it  is  assumed  that  the  error  in  Ae  is  plus  or  minus  E(ae),  the  simplest 
possible  symmetric  distribution,  then  calculated  ranges  are  the  following: 

A  A 

r  =  L  sin  e  _  R  _  R 
+  A6  +  t(\e)  7  ~  £(ao)  1  +  e 

Ad 

(47) 

o  =  L  sin  e  _  R  _  R 
-  “  A6  -  £ ( a e)  Eaq  ~  1  -  e 

1  -  A6 


where 

.  E(a«) 
e  •  -7T- 

The  range  bias  Is; 

2  A 

BIAS(R)  =  * . (48) 

1  -  e^ 

2  2 

Note  that  the  bias  is  positive  if  e  is  less  than  one,  negative  if  e  is 
greater  than  one,  and  approaches  plus  or  minus  infinity  as  e^  increases  to 
one  or  decreases  to  one,  respectively. 

The  derivations  of  the  optimum  geometric  weighting  factors  were  based 
upon  the  assumption  that  the  TMA  solution  errors  were  linear.  This  is  a  valid 
assumption  for  small  range  and  bearing  measurement  errors  at  short  ranges,  but 
it  is  not  necessarily  valid  for  the  moderate  range  and  bearing  errors  and  the 
longer  ranges  that  are  often  encountered  in  actual  target  engagements. 
Therefore,  in  order  to  determine  the  optimum  geometric  weighting  factor  under 
the  more  realistic  conditions  where  range  bias  may  be  a  significant  fraction 
of  the  total  TMA  range  error,  the  simulation  analysis  assumed  bearing  errors 
of  plus  or  minus  0.5  degrees  and  range  errors  (due  to  a  delta-tau  error  in  the 
WAA  of  plus  or  minus  2.5  microseconds)  which  resulted  in  biased  range  inputs 


to  the  TMA  algorithms.  The  results  of  this  simulation  analysis  for  the  (2, 

2  x  2;  R  -  B)  and  the  [KqPT(R  -  B)  +  (B  -  0)]  TMA  algorithms  are  given  in 
Figure  3-5,  which  shows  a(R)  and  RMS(R)  versus  for  the  two  algorithms  and 
for  each  of  the  five  run  geometries.  The  RMS(R)  is 

1 

RMS(R)  =  a2( R)  +  BIAS2 (R )  (49) 

The  optimum  geometric  weighting  factor  for  the  [KqPj(R  -  B)  +  (B  -  0)1  TMA 
algorithm  using  the  RMS(R)  curves  varies  somewhat  for  the  five  runs,  but 

is  a  reasonable  value  overall,  just  as  predicted  by  Eq.  (33).  The  optimum 
geometric  weighting  factor  for  the  (2,  2  x  2;  R  -  B)  TMA  algorithm  probably 
lies  closer  to 

3 

than  to 

4 

as  was  predicted  by  Eq.  (16).  However,  the  difference  is  so  small  that  Eq. 
(16)  may  be  used  as  the  optimum  weighting  coefficient. 

3.4  MANEUVERING  TARGET  INDUCED  TMA  ERRORS 

Maneuvering  targets  pose  a  special  problem  for  TMA  algorithms  since  they 
are  all  based  upon  Eqs.  (1)  and  (2),  which  are  for  a  constant  velocity  target 
motion.  In  order  to  directly  compare  the  accuracies  of  the  (B  -  0),  [(KqPj 
(R  -  B)  ♦  (B  -  0)]  and  (2,  2  x  2;  R  -  B)  TMA  algorithms,  a  simulation  analysis 
was  performed  using  the  nine  run  geometries  shown  in  Figure  3-6,  eight  of 
which  contain  a  target  maneuver.  Note  that  part  of  each  of  the  target 


maneuver  runs  is  the  same  as  the  constant  velocity  target  run.  Four  sets  of 
data  (at  the  beginning  and  end  of  each  leg  of  platform  motion)  were  used  to 
obtain  each  solution,  which  gave  a  total  of  256  solutions  per  run  that  were 
then  used  to  obtain  the  standard  deviation  and  bias  of  the  TMA  solution  for 
range,  down-range  target  velocity  and  cross-range  target  velocity.  This  was 
done  for  each  run  and  for  each  of  the  three  TMA  algorithms.  The  analysis 
results  are  listed  in  Table  3-1.  The  absolute  mean  values  for  the  eight 
maneuvering  target  runs  are  given  in  Figure  3-7  as  a  function  of  K  for  K  =  0, 
the  (B  -  0)  TMA  algorithm;  K  =  KqP-j- ,  the  [KgpT  (R  -  B)  +  (B  -  0)]  TMA 
algorithm;  and  K  =  1,  the  (2,  2  x  2,  R  -  8)  TMA  algorithm.  The  dashed  curves 
in  the  a(R),  BIAS(R)  and  RMS(R)  versus  K  figures  are  for  run  No.  1,  which  is 
the  non-maneuvering  target  run.  Note  that  the  range  errors  for  K  =  1,  the  (2, 
2  x  2;  R  -  B)  algorithm,  are  not  significantly  greater  for  the  maneuvering 
target  runs  than  for  run  No.  1.  Conversely,  the  RMS(R)  error  for  K  =  0,  the 
(B  -  0)  algorithm,  is  2.7  times  as  large  for  the  maneuvering  target  runs,  due 
primarily  to  the  large  BIAS(R)  error. 

While  the  eight  maneuvering  target  runs  that  were  analyzed  represent  only 
a  small  sample  out  of  the  infinite  set  of  possible  target  maneuvers,  it 
nevertheless  bridges  the  gap  between  non-maneuvering  target  runs  and  those 
maneuvering  target  runs  which  can  be  detected  by  a  maneuver  detection 
algorithm.  The  analysis  shows  that  whenever  range  measurement  data  is 
available,  the  (2,  2  x  2;  R  -  B)  TMA  algorithm  gives  the  most  accurate  range 
solution  and  a  reasonably  accurate  target  velocity  solution  for  a  maneuvering 
target.  The  (B  -  0)  TMA  algorithm  solution  for  a  maneuvering  target  tends  to 
diverge,  which  causes  unacceptably  large  bias  errors.  The  [KgpT  (r  _  b)  + 

(8  -  0)1  TMA  solution,  as  might  be  expected,  has  bias  errors  that  are 
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Run  No. 


Run  No. 


1 


0.0901 


Table  3-1.  Maneuvering  Target  Results 


(a)  K  =0,  the  (3-0)  TMA  Algorithm 

Bias  f  R )  ’<V  B(V 


B(V 


0.1065 

0.0084 

0.1865 

0.2579 

0.092 

-0.1403 

0.1225 

0.2176 

0.0939 

-0.3518 

0.0876 

-0.1007 

0.2214 

0.5823 

0.0533 

-0.196 

0.2645 

0.1909 

14 

0.266 

0.213  0.0 

16 

0.0417 

0.590  -0 .3i 

7 

1.185 

0.074  0.3 

* 

2.122 

0.356  -0.0 

18 

-0.964 

0.077  0.4i 

0.445 

0.271  0.2 

7 

3.069 

0.189  -0.2; 

!9 

0.784 

0.321  0.5- 

13 

-2.685 

0.210  -0.3' 

!9 

1.412 

0.237  0.3 

*0 

1.02 

0.161  0.1- 

-  B) 

+  (B  -  0)]  TMA  Algorithm 

Bi'as(R) 


0.0061 


2 

0.0975 

0.167 

0.681 

0.644 

0.131 

-0.118 

3 

0.0851 

-0.134 

0.681 

0.067 

0.083 

0.232 

4 

0.0891 

0.200 

0.542 

1.398 

0.131 

-0.200 

5 

0.0924 

-0.256 

0.772 

-0.995 

0.070 

0.392 

6 

0.0772 

-0.066 

0.640 

-0.438 

0.134 

0.256 

7 

0.1295 

0.307 

0.831 

1.355 

0.087 

-0.386 

8 

0.582 

-0.143 

0.467 

-0.355 

0.121 

0.182 

9 

0.1354 

0.151 

0.961 

0.076 

0.069 

-0.193 

Mean 

0.0956 

0.178 

0.698 

0.666 

0.103 

0.245 

7 

0.0241 

0.070 

0.148 

0.497 

0.027 

0.091 

Table  3-1.  Maneuvering  Target  Results  ( Continued) 


(c)  K  =  1,  the  (2,  2  x  2,  R  +  B)  TMA  Algorithm 


Run  No. 

a(R) 

nr” 

Bias(R) 

■(V 

B(V 

<l(VCR) 

8<vcr> 

ft 

VP 

VP 

1 

0.160 

0.040 

1.038 

-0.140 

0.154 

-0.046 

2 

0.107 

0.194 

0.753 

0.734 

0.137 

-0.117 

3 

0.160 

0.036 

1.052 

-0.214 

0.126 

0.132 

4 

0.156 

0.071 

0.896 

0.535 

0.177 

-0.301 

5 

0.161 

0.032 

1.075 

-0.328 

0.093 

-0.400 

6 

0.167 

0.100 

0.987 

1.039 

0.175 

0.306 

7 

0.162 

0.049 

1.059 

0.023 

0.114 

-0.453 

8 

0.148 

0.025 

0.870 

-0.289 

0.181 

0.050 

9 

0.163 

0.045 

1.082 

-0.093 

0.081 

-0.186 

Mean 

0.153 

0.069 

0.972 

0.407 

0.136 

0.243 

a 

0.018 

0.052 

0.112 

0.321 

0.037 

0.135 
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intermediate  between  the  other  two  algorithms,  but  the  standard  deviation 
errors  are  significantly  less  than  for  either  of  the  others.  Therefore,  the 
root-mean-square  errors  are  not  unreasonably  higher  than  those  of  the  (2,  2 
2;  R  -  8)  algorithm. 

The  computer  programming  and  the  computations  used  for  the  simulation 
analyses  of  the  TMA  algorithms  were  provided  by  Mr.  Juan  Arvelo  of  the  Naval 
Surface  Weapons  Center,  Code  U21.  His  significant  contribution  to  this  work 
is  gratefully  acknowledged. 


CHAPTER  4 


CONCLUSIONS  AND  RECOMMENDATIONS 

This  report  has  shown  that  the  accuracies  of  TMA  solutions  are 
significantly  improved  by  the  proper  addition  of  range  measurement  data  to  the 
least-squares  1  inear-regression  TMA  algorithms,  all  of  which  belong  to  the 
non-recursive  (or  batch)  type  of  linear  prediction  filter.  Even  though  range 
measurement  data  can  also  improve  the  accuracies  of  Kalman  filter  TMA 
algorithms  (which  are  a  recursive  type  of  linear  prediction  filter),  the 
non-recursive  type  of  filter  is  inherently  more  stable  ( i . e . ,  less  prone  to 
give  a  diverged  solution)  and  more  accurate.  In  fact,  the  recursive  filter  is 
simply  a  constant-velocity-target  tracker,  given  that  the  initial  target 
position  and/or  velocity  are  either  known  or  assumed  with  a  sufficient  degree 
of  accuracy.  Conversely,  the  non-recursive  filter  actually  solves  for  the 
constant-velocity  target  position  as  a  function  of  time  (i.e.,  the  target 
position  at  a  given  time  and' the  target  velocity).  Therefore,  it  is 
recommended  that  all  Kalman  filter  TMA  algorithms  for  use  in  underwater  fire 
control  systems  be  replaced  by  the  appropriate  least-squares  linear-regression 
TMA  algorithms. 

The  most  accurate  TMA  algorithm  for  using  both  range  and  bearing  data  was 
shown  to  be  the  [KqPj(r  -  B)  +  (B  -  0)]  TMA  algorithm.  Practical  methods  of 
determing  Kgpy  Were  derived  and  shown  to  be  valid  for  a  variety  of  typical 
target  run  geometries.  It  is  recommended  that  this  algorithm  be  implemented 
into  the  fire  control  systems  of  ships  and  submarines  that  have  accurate  range¬ 
measuring  sonar  systems. 
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Real  targets  will  most  likely  maneuver  at  frequent  pseudo-random 
intervals  when  in  potentially  hostile  areas  or  when  alerted  to  the  presence  of 
a  threat  target.  The  (2,  2  x  2;  R  -  B)  TMA  algorithm  was  shown  to  be  the 
least  affected  by  target  maneuvers.  Therefore,  it  is  recommended  that  this 
algorithm  also  be  implemented  into  the  above  fire  control  systems. 

A  (B  -  0)  TMA  algorithm  was  initially  implemented  in  the  Mk  113  Fire 
Control  System  where  it  was  called  CHURN.  It  has  been  supplanted  by  an 
operator-computer  interactive  system  called  MATE  which  is  inherently  superior 
to  the  basic  (B  -  0)  TMA  algorithm  since  it  contains  a  method  for  detecting 
target  maneuvers  and  for  casting  out  bad  data.  It  is  highly  recommended  that 
all  the  TMA  algorithms  be  provided  with  maneuver  detection  and  bad  data 
removal  algorithms.  MATE  calculates  from  the  TMA  solution  what  each  bearing 
would  have  been  at  the  time  of  each  measurement  and  subtracts  it  from  the 
measured  bearing  to  give  a  plot  of  bearing  deviation  versus  time,  which  is 
displayed  on  an  output  screen  for  the  operator  to  view.  Deviations  which 
increase  or  decrease  linearly  with  time  indicate  a  target  maneuver.  A  few 
scattered  deviations  that  are  significantly  larger  than  the  others  may  be  cast 
out  as  representing  bad  data  and  a  new  TMA  solution  can  be  calculated  with 
these  data  deleted.  An  algorithm  can  be  developed  to  automatically  cast  out 
bad  data  points  'and  relieve  the  operator  from  this  task.  It  is  also 
recommended  that  the  (8-0)  TMA  algorithm  be  implemented  for  use  whenever 
range  data  is  unavailable.  This  is  simply  done  by  setting  K  =  0  in  the 
[K(R  -  B)  +  (B  -  0)]  TMA  algorithm. 

With  only  slight  modifications  both  the  (2,  2  x  2;  R  -  B)  and  the  (B  -  0) 
TMA  algorithms  can  be  used  whenever  the  value  of  V^,  the  down-range  component 
of  target  velocity,  is  either  known  or  assumed  to  be  known  as  described  in  the 
report.  The  usefulness  of  these  algorithms  has  not  been  established. 


However,  Ekelund  ranging  is  extensively  used  in  field  tests  and  it  provides 
range  solutions  that  are  often  more  accurate  than  those  from  the  TMA 
algorithms.  Since  the  (S  -  0)  TMA  algorithm  for  an  assumed  value  of  V^, 
given  by  Eq.  (46),  is  the  least-squares  1 i near-regression  formulation  of  the 
Ekelund  algorithm,  it  should  provide  even  more  accurate  range  and  cross-range 
velocity  solutions.  A  similar  argument  holds  for  the  (2,  2  x  2;  R  -  B)  TMA 
algorithm  for  an  assumed  value  of  V^,  as  given  by  Eq.  (36).  It  is, 
therefore,  recommended  that  these  two  algorithms  be  field  tested  to  determine 
whether  or  not  they  should  be  implemented  into  the  appropriate  fire  control 
systems. 

The  optimum  weighting  coefficient,  oqPj>  for  weighting  each  linear 
equation  relative  to  the  others  was  derived  for  each  algorithm  and  then  tested 
by  simulation  analysis  for  which  the  input  data  errors  were  not  unreasonably 
large,  but  still  large  enough  to  give  biased  range  measurements.  These  biases 
affected  the  optimum  geometric  weighting  factors  somewhat,  but  not  enough  to 
justify  changing  them.  It  is  recommended  that  the  appropriate  aopT  be  used 
in  each  of  the  TMA  algorithms.  It  may  also  be  desirable  to  provide  the  option 
of  including  a  time  dependent  factor  to  agpT  in  order  to  reduce  the 
magnitudes  of  the  TMA  solution  bias  errors  that  are  caused  by  a  target 
maneuver.  An  empirically  determined  time-dependent  weighting  factor  is  the 
following: 


(50) 


where 


vp  =  platform  speed  in  knots 
fci  *  time  after  t0  in  minutes 

Q 

i  =  averaged  measured  range  at  time  ti  in  kiloyards 
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CHAPTER  5 


NOMENCLATURE 


(U,  W)  =  initial  target  position  (unknown) 

(U,  W)  =  target  velocity  (unknown) 

(U.,  W.)  *  (U  +  t.U,  W  *  t.W) 
t  =  t.  =  time  since  initial  time,  ftj  =0) 

(X,  Y)  =  (X.,  YJ  =  platform  position  at  time  t^ 

(X,  Y)  *  (Xr  Y^ )  =  measured  platform  position  at  time  t. 

E(X)  =  X  -  X  ■  error  in  measured  value  of  X 
E ( Y )  =  Y  -  Y  =  error  in  measured  value  of  Y 
R  =  R^  =  target  range  (from  platform)  at  time  t^ 

R  =  R^  *  measured  target  range  at  time  t. 

E(R)  ■  R  -  R  *  error  in  measured  range 

B  =  B.  =  true  bearing  to  target  at  time  tj 
B  *  B.,-  3  measured  true  bearing  to  target  at  time  t-j 
E(B)  -  B  -  B  =  error  in  measured  true  bearing 

s  =  sin  B  =  sin  [B  -  E(B)]  a  sin  B  -  (cos  B)  E(B)  ~  S  -  cE(B) 

s  =  sin  B 

c  =  cos  8  =  cos  [B  -  E(B)]  »  cos  B  +  (sin  B)  E(B)  =  c  +  s  E(B) 
c  =  cos  B 

E ( s )  a  (cos  B)  E(B)  =  c  E(B)  =  error  in  sin  B  due  to  E(B)  (a  first  order 
approximation) 

E ( c )  «  -  (sin  B)  E(B)  =  -  s  E(B)  =  error  in  cos  B  due  to  E(B)  (a  first 
order  approximation) 


a  =  aj  =  a  weighting  coefficient  for  the  ith  equation  in  the  set  of 
n  equations 

K  =  a  weighting  coefficient  for  the  range-bearing  TMA  solution  relative 
to  the  bearings-only  TMA  solution 

VCR  =  VCR  ^  =  t^)  -  cross-range  component  of  target  velocity  at  time  tk 

VR  -  VR  (t  -  tk)  =  down-range  component  of  target  velocity  at  time  t^ 

»  (Vp  cos  e)^  *  down-range  component  of  platform  velocity  at  time  t 
K 


PLATFORM, 


RpiATFORM  =  mean  down“ran9e  component  of  platform  velocity  from  time 
t  =  0  to  time  t  *  t|< 

(R)k  =  TMA  range  at  time  t^ 

MEDIAN(R)  =  median  value  of  measured  range 
A, 

A 

A 
A 
A 
A 
A 


'll = 

T.c2 

B11  = 

[\2  * 

Z  ate2 

B^2  *  Zat 

L13  = 

Z  aSC 

B12  =  Z  at2 

l14  = 

Z  atSC 

l22  = 

Zat2C2 

F1  =  Z  o(Rc  +  X) 

‘24  = 

Z  at2SC 

F2  *  ZatfRe  +  X) 

‘33  = 

ZaS2 

F  3  =  ZafRs  +  Y) 

‘34  = 

ZatS2 

F4  =  Z  at(R$  +  Y) 

‘44  = 

Zat2S2 

'.■I 

,aC(R  +  CX  +  sY) 

Ej  =  Z  ats(  sX  -  cY) 

'2-1 

-  atc(R  +  cX  +  sY) 

E2  =  Zats(sX  -  cY) 

'3=1 

,as(R  +  cX  +  sY) 

E3  =  - Z  aC ( sX  -  cY) 

°4  =  1  ats(R  +  CX  +  sY) 


E4  =  -  _  atx ( sX  -  cY) 


APPENDIX  A 


Kqpt  For  [K(R  -  B)  +  (8  -  0)1  TMA  Algorithm  Using  D/E  or 
Multipath  Ranging  Data 


Figure  A-l  shows  the  geometry  for  D/E  or  multipath  ranging.  The  figure 
shows  straightline  propagation  paths,  when  in  reality  these  paths  may  be 
curved  due  to  temperature  and  pressure  gradients.  Therefore,  it  is  assumed 
that  each  actual  path  has  been  ray-path-curvature  corrected  to  its  equivalent 
straightline  path.  The  target  range  is 

R  =  (D  -  P)  cot  4  +  (D  -  T)  cot  (4  +  2a)  (A-l) 

where 

D  =  ocean  depth  at  the  reflection  point 
a  =  ocean  bottom  down-range  slope  at  the  reflection  point 
=  E( a)  if  a  is  not  known 
P  =  platform  depth 

4  -  vertical  arrival  angle  of  the  broadband  radiated  target  signal 
T  =  target  depth 

The  range  error,  to  a  first  order  approximation  is 

E(R)  =  [cot  4  +  cot  (4  +  2a)]  E(D)  -  cot  (4  +  2a)  E(T) 


E(d)  +  24D  -  E(a) 

sin  ( 4  +  2a) 


(A-2) 


If  it  is  assumed  that  a  <<  4  and  P,  T  <<  D 


r,n,  .  „  2  E(D)  -  E(T)  ,  2  E ( a)  -  2  E(«5) 

E(R)  r  — - +  — nrra — 


a(R)  =  R  Va2  (R)d  +  a2  (R)y  +  a2  (R)  a  +  J  (R)  ^ 


(A-3) 


where 


r 

17' 


tform 


Surface 


T 

K 
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Refl  ected 
Surface 


Figure  A-1.  Geometry  for  D/E  or  Multipath  Ranging. 
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a(T) 


*(R)a  =  sTTT?  a(a)  = 


R  ct( a) 

sin  t  cos  1 6 


o(R) j  = 


2R 


i  ~  si  n  2i 


:d)  = 


R  a  ( a) 

sin  6  cos  t 


The  value  of  K  which  minimizes  the  [K(R  -  B)  +  (3  -  0)]  TMA  solution  errors  is 


K0PT 


( A— 4 ) 


A  reasonable  example  of  the  errors  a(8)  and  those  which  comprise  a(R)  are 
the  following: 

Let 


o(D)  =  450'  =  0.15  kyds 
a(d)  *  1°  =  0.0175  radians 

cr(a)  -  1° 

o(B)  *  a ( a )  *  1°  (in  bottom-bounce  mode) 

Substituting  these  standard  deviation  errors  into  Eqs.  (A-3)  and  (A-4)  gives 


q(R) 

R 


(0.158) 


0.01236) 


/ 


kopt 


3.056  10'4  D2 


0.025  +1.53  10"4 


+  4  D2\2 

*  V 


(A-5) 


Figure  A-2  shows  as  a  function  of  range  for  four  different  ocean  bottom 
depths.  This  set  of  curves  allows  an  operator  to  easily  obtain  a  reasonable 
value  for  Kgpy,  given  that  the  assumed  errors  are  reasonably  accurate.  If 
other  values  of  the  assumed  errors  are  known  to  be  more  accurate,  they  can  be 
substituted  into  Eqs.  (A-3)  and  (A-4)  and  another  set  of  Kqpj  curves  as  a 
function  of  range  and  ocean  depth  can  be  plotted  for  use  by  the  TMA  operator. 
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APPENDIX  B 


LEAD-LAG  VERSUS  LAG-LEAD  PLATFORM  TACTICS 

Bearings-only,  (B  -  0),  TMA  solution  errors  are  known  to  be  dependent 

/ 

upon  the  geometry  of  the  target  and  platform  motion.  The  (B  -  0)  TMA 
algorithm  is  based  upon  the  assumptions:  (1)  that  the  target  speed  and  heading 
(i.e.,  velocity)  are  constant  throughout  the  solution  time  interval  and  (2) 
that  the  platform  maneuvers  (i.e.,  changes  speed,  heading  or  both)  at  least 
once  during  the  solution  time  interval.  It  has  been  shown  that  the  most 
accurate  TMA  solution  for  a  given  platform  speed  and  solution  time  interval  is 
obtained  when  the  platform  conducts  a  lead-leg  followed  by  a  turn  and  then  a 
lag-leg  of  motion  (i.e.,  a  lead-lag  tactic).^  A  lead-leg  of  platform  motion 
is  defined  as  one  in  which  the  platform  moves  in  the  same  cross-range 
direction  as  the  target,  and  a  lag-leg  is  one  in  which  the  platform  moves  in 
the  opposite  cross-range  direction. 

An  example  will  be  used  to  show  why  the  lead-lag  tactic  gives  more 
accurate  range  estimates  than  the  lag-lead  tactic.  Figure  B-l  shows  three  TMA 
run  geometries:  one  in  which  the  target  is  moving  cross-range  at  the  same 
speed  as  the  platform  and  the  platform  is  conducting  a  lead-lag  tactic; 
another  in  which  the  target  motion  is  the  same  as  the  first  run,  but  the 
p^tform  is  conducting  a  lag-lead  tactic;  and  one  in  which  the  target  is 
stationary  and  the  platform  conducts  the  same  two-leg  tactic  as  either  of  the 
first  two  run  geometries  (since,  by  symmetry,  they  are  equivalent  for  a 
stationary  target).  The  target  range  is  about  twelve  times  the  length  of  each 
leg  of  platform  motion  for  all  three  run  geometries  (i.e.,  R  »  L  *  VpT). 


The  TMA  algorithm  solves  for  the  X  and  Y  components  of  target  velocity 
and  initial  position.  These  Cartesian  coordinate  solutions  can  be  used  to 
calculate  target  range  at  any  time  before,  during  or  after  the  solution  time 
interval  (which  for  these  examples  is  from  t  =  0  to  t  =  3T,  where  T  is  the 
duration  of  each  leg  of  platform  motion  and  of  the  turn).  The  TMA  range  error 
is  the  difference  between  this  calculated  range  and  the  actual  target  range. 
The  same  analytic  method  as  was  previously  described  was  used  to  obtain  the 
standard  deviation  range  error,  a(R),  as  a  function  of  a  given  standard 
deviation  bearing  measurement  error  for  the  three  TMA  run  geometries  shown  in 
Figure  B-l.  The  results  are  shown  in  Figure  B-2  which  shows  relative  standard 
deviation  range  error  as  a  function  of  time.  The  a(R)  is  normalized  relative 
to  the  minimum  o(R),  which  is  the  same  for  all  three  run  geometries  but  which 
occurs  at  different  times.  The  minimum  crC R)  for  the  stationary  target  occurs 
at  the  midpoint  (t  =  1.5T)  of  the  solution  time  interval,  and  the  a(R)  at  the 
beginning  and  at  the  end  of  the  solution  time  interval  is  twice  as  large. 

Note  that  the  a(R)  curves  for  the  three  run  geometries  are  identical  except 
for  a  time  shift.  The  a(R)  curve  for  the  cross-range  moving  target  is  shifted 
in  time  by  minus  2T  for  the  lag-lead  platform  motion  tactic  and  by  plus  2T  for 
the  lead-lag  platform  motion  tactic.  This  has  the  effect  of  making  the 
lead-lag  TMA  range  solution  more  accurate  at  a  future  time,  3.5T,  (or  0.5T 
after  the  end,  3T,  of  the  solution  time  interval!  than  at  any  time  within  the 
solution  time  interval.  It  also  has  the  effect  of  making  the  lag-lead  TMA 
solution  range  error  at  time  0.5T  after  the  end  of  the  solution  time  interval 
4.6  times  as  large  as  that  of  the  lead-lag  TMA  solution.  The  range  error  of 
the  lead-lag  TMA  solution  at  the  end  of  the  solution  time  interval  (t  »  3T)  is 
about  0.57  that  for  the  stationary  target  TMA  solution  and  about  0.28  that  for 
the  lag-lead  TMA  solution. 
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Figure  B-l .  Platform  Target  Encounter  Tacti 
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If  the  target  cross-range  velocity  had  been  greater  (or  less)  than  the 
platform  cross-range  velocity,  the  time  displacements  of  the  a(R)  curves  would 
also  have  been  greater  (or  less).  This  example  illustrates  very  succinctly 
the  superiority  of  the  lead- lag  over  the  lag- lead  platform  motion  tactic  for 
2  bearings-only  TMA  algorithms,  since  the  difference  in  their  accuracies  is 


a  function  of  the  run  geometry. 
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